Magnetoresistance in bilayer graphene via ferromagnet proximity 

effects 



Y. G. Semenov, J. M. Zavada, and K. W. Kim 
Department of Electrical and Computer Engineering, 
North Carolina State University, Raleigh, NC 27695-7911 

Abstract 

A drastic modification of electronic band structure leading to the magnetoresistance is predicted 
in bilayer graphene when it is placed between two ferromagnetic insulators. Due to the exchange in- 
teraction with the proximate ferromagnet, the electronic energy dispersion in the graphene channel 
strongly depends on the magnetization orientation of two ferromagnetic layers, Mi and M2. While 
the parallel configuration Mi = M2 leads to simple spin splitting of both conduction and valence 
bands, an energy gap is induced as soon as the angle 9 between Mi and M2 becomes non-zero with 
the maximum achieved at = vr (i.e., antiparallel alignment). Consequently, bilayer graphene may 
exhibit a sizable magnetoresistance effect in the current-in-plane configuration. A rough estimate 
suggests that the resistance change on the order of tens of percent is possible at room temperature. 
This effect is expected to become more pronounced as the temperatures decreases. 
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With the advent of free-standing atomically thin graphite films (or graphene),- unusual 
properties related to the two-dimensional Dirac-like relativistic spectrum of the honeycomb 
carbon lattice^ put graphene in the forefront of emerging carbon based electronics. Among 
others, the half-integer quantum Hall effect^i^i^ observed even at room temperature,^ high 
carrier mobility, easy control of electron and hole concentrations via variation of applied 
bias, absence of weak localization and universal minimal conductivity at the Dirac point 
with zero density of states have attracted significant theoretical and experimental attention 
to electronic transport in graphene .-i^>^ 

The spin dependent properties of graphene also offer fascinating opportunities. In most 
studies, the starting point of consideration is its extremely small spin-orbital coupling com- 
pared to typical semiconductors.— Consequently, graphene exhibits long electron spin- 
relaxation time^ and mean free pathsF^^ even at room temperature.^^ However, this very 
advantage (i.e., the weak spin-orbital interaction) presents a severe challenge in spin ma- 
nipulation or selection via electrical control. One possible approach may be to utilize the 
electron exchange interaction with proximate ferromagnetic layers and the resulting effective 
magnetic field.— 

It was found very recently that bilayer graphene (BLG) offers a host of novel phenom- 
ena through external modification of the energy band structure.— ii^i22i2ii2^ Particularly, it 
was shown that a gap opens up between the conduction and valence energy bands when 
a potential difference u is introduced in the two graphene layers. Moreover, the parabolic 
band structure near the non-equivalent K and K' points^ transforms to a Mexican-hat-like 
dispersion with u 7^ 0.— Likewise, a non-trivial response may be expected in BLG to the 
symmetry-breaking spin- dependent interactions. 

In this study, the properties of BLG sandwiched between two insulating/ dielectric ferro- 
magnets with magnetic moments Mi and M2 are exploited based on an 8 x 8 tight-binding 
model. Similar to the electrical bias, the calculation predicts a significant modification of 
the electronic band structure when the constituting graphene layers are subject to un-equal 
exchange interactions with the proximate magnetic ions. A particularly interesting fea- 
ture appears for ferromagnetic layers with identical magnitude of the magnetic moments 
(Ml = M2), where u = and the potential asymmetry is not a factor. As the alignment of 
Ml and M2 deviates from each other, the spin interactions reveal an asymmetry that can 
alter the BLG energy bands (including an energy gap) and subsequently the in-plane con- 
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ductivity. The resulting magnetoresistance effect can play an important role in the carbon 
based spintronics. In the case of monolayer graphene, on the other hand, the net effect of 
two magnetic layers simply reduces to electron spin splitting in the effective field induced 
by the vectorial sum Mi + M2. 

Figure 1 schematically illustrates the specific structure under consideration. It resembles 
the ferromagnet-metal hybrid structures'*^ that reveal a giant magnetoresistance owing to 
the spin-dependent conductivity. The bottom ferromagnetic dielectric layer (FDL) possesses 
the magnetization Mi that can be pinned along the direction of the x axis by an antifer- 
romagnetic substrate. The top FDL may be constructed from the same material but its 
magnetization vector M2 can be rotated on the x-y plane (by an external magnetic field) 
forming an angle 6 with Mi. The influence of FDL magnetization on BLG electronic struc- 
ture can be realized in actual structures through either the direct exchange interaction with 
magnetic ions (assuming an overlap between the carbon vr-orbitals and unfilled shells of 
the magnetic ions in FDLs) or an indirect interaction via the ligands of FDLs. Thus, the 
problem can be modeled in the mean field approximation with the Hamiltonian 

H = Hbl + Pi«MiS+P2«M2S , (1) 

where Hbl is the spin-independent BLG Hamiltonian. Two remaining terms of Eq. ([T]) 
describe the energy of an electron spin S in the effective fields (in units of energy) aMi 
and aM2 of the proximate FDLs. Accordingly, projection operator Vi {V2) is 1 for the 
electron localized at the bottom (top) carbon monolayer and otherwise. Parameter a is 
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proportional to the carrier-ion exchange constant as evaluated in Refs. 

In the case of low energy electronic excitations, the tight-binding approximation-^ can 
accurately describe the BLG band spectra as recently demonstrated in a density functional 
calculation.— Hence, we adopt the tight-binding Hamiltonian near the valley extrema in a 
basis that constitutes the components {Ai |, Ai |, B2 t? -B2 i, ^2 T) ^2 I, Bi t, Bi |) for 
the K valley and {B2 T, B2 i, Ai T, ^1 i, t, B^ i, A2 T, ^2 i) for the K' valley. Here, 
Ai and B^ correspond to the electron amplitudes at inequivalent sites of the bottom {i = 1) 
and top {i = 2) graphene layers as shown in Fig. 1(b), and | and | denote spin up and spin 
down states. 

The Hamiltonian Hbl includes the lateral coupling for nearest carbon atoms in the bot- 
tom [Ai — Bi) and top {A2 — B2) layers with the matrix element 7 (= 3 eV) as well as the 
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interlayer coupling for A2 — Bi and Ai — B2 dimers with the matrix elements 71 (= 0.4 eV) 
and 73 (= 0.3 eV), respectively; hence, 73 < 71 -C 7. At zero magnetic field, Eq. ([1]) for the 
lowest electronic states in the K and K' valleys can be expressed in terms of the in-plane 
electron momentum k = [k^, ky) defined from the centrum of each valley as 
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the Fermi energy in monolayer graphene.- In addition, the term ^3 = v/3a73/2^ v) 
is responsible for the trigonal warping, where a = 0.249 nm is the length of lattice unit 
vector. The secular equation for Eq. ([2]) is solved under the conditions of (i) zero bias and 
(ii) identical top and bottom FM materials (i.e., u = 0, Mi = M2 = M). For simplicity, 
the transfer matrix elements for off-center sites Ai and B2 are also ignored {v^ 0). This 
approximation does not change qualitatively the results for the energies E > 1 meV as 



discussed in Ref. 
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The specified conditions leads to the energy spectra of each valley consisting of eight 
non-degenerate branches en{k) that are identical for conduction and valence bands and 
isotropic with respect to the valley centrum (i.e., independent of Two spin pair solutions 
correspond to the excited states with energies |en(^)| ^ 7i that are beyond the current 
interest. Hence, only the remaining four low-energy bands e„ = en{k) are considered. As 
it is convenient to normalize the parameters in units of 71, the dimensionless momentum 
p = vk/'ji and the exchange field G = q;M/7i are introduced hereinafter. The low energy 
bands = ±'yiE± can be expressed as 
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for ^ 6' ^ vr. Two solutions with e„ > correspond to the lowest conduction bands, while 
their mirror images with respect to the zero energy describe the highest valence bands. 

The most remarkable outcome of the calculation is the presence of an energy gap Eg 
between the lowest conduction band and the highest valence band for 7^ as shown in 
Fig. 2. When the orientation of Mi and M2 is in parallel alignment [9 = 0, Fig. 2(a)], 
the net effect of the exchange interaction simply lifts the two-fold spin degeneracy resulting 
in two pairs of spin-split bands that cross each other at p = ^/G{1 + G/2)/2. However, 
once they are disaligned, the electronic bands become of mixed spin character (e.g., with 
both parallel and antiparallel components to the x direction). Subsequent anti-crossing 
opens up the gap that progressively grows with 9. At 9 = n [i.e.. Mi = — M2, Fig. 2(d)], 
the conduction and valence bands are merged to form two doubly degenerate states with 
the maximal Eg = G/y/l + CP a.t p = Gy2 + CP/2\^l + CP. This spin degeneracy is not 
surprising because the equivalence of top and bottom graphene layers makes the Mi and 
M2 (= —Ml) directions indistinguishable [i.e., the dependence of Eqs. and (jl]) on sign 
in ±G disappears]. Note that in the case of monolayer graphene, the effect of two FDLs 
in a similar configuration M2 = —Mi gives rise to net cancellation of spin-dependent band 
modification. 

In general, the bandgap induced in BLG can be expressed in term of the deviation angle 
9 as 



Evident from this equation, the strength of the effective field G determines the size of Eg. 
If its magnitude is comparable to the thermal energy (strictly speaking 71 G ^ ksT as G is 
a normalized quantity), the orientation dependence in Eq. ([5]) may manifest itself through 
the variation of electron/hole population. At the same time, the carrier velocity in BLG 
is also strongly affected as indicated by the flattening of the bands with increasing 9 in 
Fig. 2. Clearly the conductivity in the BLG channel can be modulated by the effective 
fields from the two ferromagnetic barriers, leading to a sizable magnetoresistance effect in 
the current-in-plane configuration. 

The effect of both conduct ion/ valence band separation and deformation on the magne- 
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toresistance R{0) can be taken into account in the one-electron approximation in terms of 
Kubo-Greenwood formula for conduct ivity^i^^ 



m,n k,k' ^ ^ '".k 
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where e is the electron charge, m and n the subband indices, Vx = dx/dt = i[H^,x]/h 
the velocity operator, and / the Fermi-Dirac function. The case of ballistic conductivity 
is examined at finite temperature.— The specific property of interest is the relative change 
o'xx{(^)/o'xx{0) as a function of 6. Utilizing G ^ 1 in most cases, our analysis reveals that it 
depends most sensitively on a single parameter 'yiG/ksT as expected. Figure 3 provides the 
calculated cr^x(6')/cr^x(0) as well as the corresponding change in the electron/hole concentra- 
tion n{9)/n{0) for two different values of '-fiG/ksT. The comparatively weaker decrease of 
n (in reference to axx) as 6 rotates signifies a major contribution of the mobility variation 
due to the flattened band structure. 

Assuming a sensitive response of M2 to the external magnetic field, one can de- 
fine the magnetoresistance as a difference of the BLG resistance -R(O) ~ l/axx{0) at 
Ml = M2 and R{tt) ~ l/axxi'n') at Mi = — M2. Figure 4 presents the normalized re- 
sult ^ = [-R(vr) — _R(0)]/-R(0) as a function of ■yiG/kBT. Since this quantity depends on 
the position of the electro-chemical potential fi as well, four different values of /ksT are 
considered. All the cases show essentially the same characteristics with the largest effect 
exhibited at yU = (i.e., fi at the mid gap). ^ also saturates at smaller effective fields as 
gets closer to the conduct ion/ valence band edge. An important point to note is that the 
magnetoresistance effect of tens of percent is possible once the exchange field strength 71G 
(= aM) is comparable to the thermal energy. By comparison, the alternative approach 
based on spin-valve devices in monolayer graphene reveals only a feeble magnetoresistance 
due to the weak dependence of the graphene conductivity on the electronic details of the 
ferromagnetic leads.— 

As a practical matter, the strength of the exchange field G is of major importance. 
Considering a large variety of potential FDL/graphene interfaces and the current lack of 
relevant information, however, a priori evaluation of G may be possible only in a very rough 
manner. Following the approach presented in Ref. Il7|, we estimate the exchange energy 
for a graphene electron interacting with the nearest stratum of magnetic ions as aM = 
{n2FM/n2c)SMJ, where n2c and n2FM are the areal concentrations of graphene carbon atoms 



and magnetic ions in the FDL, respectively, Sm is the mean value of the magnetic ion spins, 
and J is the exchange constant. The latter was found to be J = 15 meV as deduced from an 
experiment in an EuO/Al structure with Curie temperature = 69 K.^^ This constant J 
may actually be larger in the case of high temperature ferromagnets with stronger spin-spin 
inter-ions interaction. With the provision that J ~ Tc, n2FM/n2c = 0.1 — 0.2, Sm = 1 — 2.5 
and Tc = 500 — 600 K, we find 'jiG (= aM) ~ 15 — 65 meV. Although the accuracy of 
this estimate is limited, it nonetheless shows the possibility of prominent magnetoresistance 
effect even at room temperature. Lower temperatures will make its manifestation far more 
apparent for easier detection. 

It is noted that the idealized model of Eq. ([T]) may need to be expanded for detailed 
analysis of a specific structure. Imperfect interfaces between the FDL and BLG can be 
a source of a random potential that affects the electronic band structure and transport 
properties of graphene.— In addition, any difference between the top and bottom FDLs 
results in a non-zero bandgap even at parallel orientation of their magnetic moments. At the 
same time, the non-trivial manifestation of electron-electron interaction^'^ may interplay 
with the studied effects at sufficiently low temperatures. These factors (and many others)^ 
may need to be taken into consideration as the situation requires. 

This work was supported in part by the US Army Research Office and the FCRP Center 
on Functional Engineered Nano Architectonics (FENA). 
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FIG. 1: (a) Schematic illustration of bilayer graphene (two closely set planes) sandwiched between 
ferromagnetic dielectric layers of magnetization Mi and M2 (separated by the angle 6). The 
reference frame is chosen so that Mi is along the x axis; the grey arrow shows the probing current 
Ix through the graphene channel, (b) Fragment of bilayer graphene as two hexagons (view from 
above) with carbon atoms located at the vertices. Lattice sites Ai and Bi (dashed hexagon) refer 
to the bottom layer, while A2 and B2 (solid hexagon) the top layer. 7, 71, and 73 represent the 
matrix elements for electron transfer between the nearest in-plane, vertical, and slanted inter-layer 
carbon atoms, respectively. 
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FIG. 2: Electron energy spectra of bilayer graphene near the K point at different angles 9 between 
magnetization vectors Mi and M2. The reference for energy E = corresponds to the center of 
the bandgap [(b)-(d)] or the point of contact between the conduction and valence bands in the case 
^ = [(a)]. The bands are doubly degenerate aX 9 = ir [(d)]. 
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FIG. 3: Conductivity a (solid lines) and carrier concentration n (dashed lines) vs. 9 in bilayer 
graphene at two different ratios of the thermal energy fc^T to the exchange field jiG (i.e., 
ksT/^iG = 0.6 or 0.4). The electro-chemical potential is fixed at zero (i.e., at the middle of 
the energy gap). 
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FIG. 4: Magnetoresistance $, in the anti-parallel alignment (Mi = — M2) as a function of ^iG/ksT 
calculated at different locations of the electro-chemical potential ^i. The normalized value ji/kBT 
is provided for each curve. 
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